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By using Monte-Carlo implementations of fcT-factorization formula with running- 
coupling BK unintegrated gluon distributions for nucleus-nucleus collisions, we com- 
pute higher order harmonic moments of the initial density distribution for both 
RHIC(Au-HAue200GeV) and LHC(Pb-fPb@2.76TeV) collisions. We study their sensitiv- 
ity to the size of the valence parton distribution in the nucleon. 

I §1. Introduction 

> . 

Q , Spacetime evolution of the hot and dense system created from heavy ion col- 

' hsions at Relativistic Heavy Ion Colhder (RHIC) or Large Hadron Colhder (LHC) 

may be followed by hydrodynamical simulations. Indeed, the nearly perfect fluid 
picture appears to explain large elliptic flow discovered at RHIC and LHC. Recently, 
it was reported the effects of initial condition for hydrodynamics on the higher order 
flow harmonics in Au-|-Au collisions at RHIC-L' by comparing various hydrodynam- 
ical models with different initial conditions.'^'® Combined analysis of higher order 
^ ' flow harmonics may provide strong constraints on both the properties of QGP and 

Ch^. the initial higher order moments of the density distribution of gluons in coordinate 

spacc;^' It was pointed out in Ref.l^ that all moments have the same magnitude 
' when there is no length scale for nucleons inside a nucleus, while their magnitude 

decreases with the harmonic number if a length scale is introduced. In this work, 
I we compute higher order harmonic moments within Monte-Carlo version of kt fac- 

OO ■ torization formulation with running coupling BK equation, and study the effect of 

Gaussian width on the harmonic moments. 
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§2. Theoretical Model 

We will use the Monte-Carlo implementation of kt formula with unintegrated 
gluon distribution function from numerical solutions of running coupling Balitsky- 
Kovchegov (MCrcBK^® equation for the computation of gluon production in 
5^ I heavy ion collisions. This model is a extension of the Monte-Carlo KLN (MC-KLN) 

model.® First, the nucleons in the two incident nuclei are randomly sampled ac- 
cording to the Woods-Saxon distribution which simulates the effects of fluctuations 
of position of hard color sources. At each grid point, we compute the thickness func- 
tion TA(r_L) to obtain the local saturation scale and then compute the gluon density. 
Within a Gaussian nucleon approximation, thickness function (for the large-x valence 
partons) is given by 

Tp{r) = ^eM-rVm] ■ (2-1) 



typeset using VTpT^.ds (Ver.0.9) 
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Note that the mean square radius of the valence parton distribution corresponds to 
(r^) = 2B. The probabihty of nucleon-nucleon colhsion P{b) at impact parameter b 
is 

P{b) = 1 - exp[-kTpp{b)], Tpp{b) = J dhTpis) Tp{s - b) . (2-2) 

where (perturbatively) k corresponds to the product of gluon-gluon cross section and 
gluon density squared. We fix k so that integral over impact parameter becomes the 
nucleon-nucleon inelastic cross section a at at at the given energy: 

cJnn{^s) = j dH (1 - exp[-A:( V^) Tpp{b)]) , (2-3) 

with aNN = 61.36 mb for yJSNN = 2.76 TeV. Hence, P{b) broadens with increasing 
energy, even as the size B of the hard valence partons is fixed. 

In this work, we vary the width of the Gaussian within the range B = 0.2, 0.3, 0.4 
fm^. Our model applies the fc^-factorized formula in the transverse plane perpendic- 
ular to the beam axis locally. The number distribution of produced gluons is given 
by 

dNg N, f d^Pl_ 



d'^r±dy — 1 7 



d k± as (Pa{xi, +P±) ) (l)Bix2, ik± - p±) ) . 



(2-4) 

In MCrcBK,l3 ^ is obtained from the Fourier transform of the numerical results of 
the running coupling BK (rcBK) evolution equation.^i^ Note that in ()2-4p the small- 
X evolution is not treated stochastically and that the density of produced gluons 
fluctuates only due to fluctuations of the large-x valence charges. 
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§3. Initial Harmonic moments 

In Fig. [H we plot the event average of eccentricities En {n = 2, 3, 4, 5) defined 



(r^ cos(n(/)))2 + (r^ sin(n(/>))2 

where r'^ = +1]"^ ^ x = r cos (p, y = r sin cp. {■ ■ ■) means the average in the transverse 
plane. Since the eccentricities may be proportional to the magnitude of flows, it is 
important to know the initial value of these quantities. The result from Monte-Carlo 
Glauber model in which eccentricity is calculated based on the positions of point-like 
participant nucleons that is labeled by Glauber (point-like) is shown in full squares. 
Monte-Carlo Glauber results assuming hard disc nucleon are plotted in open squares. 
In the hard disc approximation, the thickness function at each grid is obtained by 
counting the number of participant nucleons as TA{r±) = number of nucleon within 5 ^ 
where the smearing area 5" = 42 mb, at top RHIC energy. The finite size of the nu- 
cleons does not affect £2 significantly, except for very peripheral collisions. However, 
Glauber with point-like nucleons yields larger values for higher moments (n > 3) 
which is consistent with the results in Ref.^^ 



Harmonic moments 
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Fig. 1. Comparison of the centrality dependence of higher order harmonic moments for Au+Au 
collisions at ^sjvjv ~ 200 GeV. 



On the other hand, one observes that the third and fifth harmonics are similar 
between the MC-GIauber (hard disc) and the MCrcBK models except for small value 
of Gaussian width B. in our model as pointed out in Refs.'^^J'El' 

We should mention that the prediction of £3 from a newly developed event gen- 
erator DIPSY based on a dipole model is larger than that of the MC-KLN modelj^ 
while DIPSY prediction for £2 is the same as MC-KLN result. This may due to the 
additional fluctuations from BFKL cascade in DIPSY which is not included in our 
models in this work. 

Figs. [2] show the harmonic moments for Au+Au collisions at ^/snn = 200 GeV 
and Pb+Pb collisions at ^/sNN = 2.76 TeV from the MCrcBK model. One ob- 
serves that as Gaussian width decreases, higher order harmonic moments increase. 
Therefore, it is important to check the sensitivity of higher flow harmonics to the 
length scale introduced by valence parton distribution in order to extract detailed 
information on the properties of quark-gluon plasma. 

§4. Summary 

We have presented results for higher order harmonic moments from the Monte 
Carlo version of kt factorization formula with rcBK small-x evolution (MCrcBK). 
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Fig. 2. Harmonic moments for central Au+Au collision (6 < 3 fm) at -^sjvjv = 200 GeV (left) and 
Pb+Pb collision at ysjviv = 2.76 TeV (right) from MCrcBK with various Gaussian width. 



The present simulations account only for the fluctuations of the valence partons in 
the transverse plane. We have shown the length scale dependence of the harmonic 
moments of the initial density distribution. It will be interesting to see how this 
affects higher order hydrodynamical flows. Also, in the future it would be interest- 
ing to systematically explore the effect of additional fluctuations from the BFKL 
evolution ladders.'^ 
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